METHOD AND SYSTEM FOR SEARCHING FOR RELATIONSHIPS BETWEEN 
BASE SEQUENCES IN GENES 



BACKGROUND OF THE INVENTION 

1 -Field of the Invention 

The present invention relates to a method and a system for searching for 
relationships between base sequences in genes. More particularly, the present invention 
relates to a method and a system for searching for the relationships between a sample base 
sequence and known base sequences based on restriction fragment length polymorphisms 
(RFLP). 

2. Description of Related Art 

An initial bacterial classification system has been constructed based on bacterial 
cell morphology, physiological properties, biochemical activities and the like. For the 
purpose of identifying an unknown bacterium based on this initial bacterial classification 
system, a plurality of physiological test results for a bacterium being identified are 
automatically recorded. The recorded physiological test results are compared with the 
physiological test results of known bacterial strains already recorded in a stored database. 
A search is made for the relationship between the unknown bacterium and known bacterial 
strains to infer the species name of the unknown bacterium. The identification method 
described above has been utilized in medical and clinical fields as well as in a wide variety 
of industrial fields such as cosmetics, food hygiene, quality control, and environmental 
hygiene. 

However, precise definitions of species for eubacteria and archaeobacteria, which 
propagate by cell division, are difflcuh to construct, in contrast to eucaryotes, which 
reproduce sexually, and for which precise definitions have been established. Therefore, 
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recently, a new classification system based on ribosomal RNA as a molecular clock has 
been constructed. To be more precise, in the new classification system, the bacteria are 
classified in accordance with differences in bacterial 16S rDNA sequences. As a result, 
there are large discrepancies between the results by the identification method described 
above and those according to 16S rDNA sequences. 

To determine the base sequences, a method utilizing a sequencer device has been 
widely used. However, the sequencer device is expensive and the number of bases that 
one sequencer device can read is limited. Thus, the method utilizing a sequencer device 
is not efficient and incurs high costs for determining the base sequences. 

Alternatively, a method has been used which involves digesting an unknown 
DNA molecule by restriction enzymes and inferring the base sequence thereof based on 
the lengths of resulting fi-agments, so as to determine the base sequence. However, this 
method requires a great amoimt of time and labor to infer the base sequence by restriction 
pattern matching of the unknown DNA molecule and known DNA molecules. Thus, this 
method is also inefficient for identifying base sequences. 

SUMMARY OF THE INVENTION 
The present invention has been made in view of the above problems, and an 
object thereof is to provide a method and a system for searching for the relationships in 
base sequences in order to efficiently analyze the relationships of an unknown base 
sequence to known base sequences based on restriction fragment length polymorphisms 
(RFLP) thereof. 

To solve the above problem, a system for searching for relationships between 
base sequences in genes according to the present invention comprises a theoretical value 
calculating portion which calculates theoretical restriction fragment length values of 
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known genes digested by restriction enzymes based on restriction enzyme data in which 
kinds of restriction enzymes, recognition sites recognized by the restriction enzymes, and 
restriction sites at which the restriction enzymes cut in the recognition sites are stored in a 
related manner, and base sequence data in which kinds of known genes and base 
sequences of the known genes are stored in a related manner, and which outputs 
theoretical restriction fragment patterns in which the kinds of known genes, the kinds of 
restriction enzymes, and theoretical restriction fragment length values are stored in a 
related manner; and a comparing portion which compares the theoretical restriction 
fragment patterns with measured restriction fragment patterns in which the kinds of 
restriction enzymes and measured restriction fragment length values obtained as a result of 
digesting a sample by the restriction enzymes followed by measurement are stored in a 
related manner, calculates the degree of similarity of the known gene to the sample, and 
outputs the similarity as analysis result data. 

The system for searching for relationships between base sequences according to 
the present invention may fiirther comprise a displaying portion that diagrammatically 
displays samples having a high degree of similarity therebetween or a knovm gene and a 
sample having a high degree of similarity therebetween in a related manner based on the 
analysis result data. 

In the system for searching for relationships between base sequences in genes 
according to the present invention, the displaying portion may display samples having a 
high degree of similarity therebetween or a known gene and a sample having a high 
degree of similarity therebetween in a related manner by a dendrogram. 

In the system for searching for relationships in base sequences according to the 
present invention, the comparing portion may calculate the degree of similarity based on 
the theoretical restriction fragment patterns and the measured restriction fragment patterns 
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by using an unweighted-pair-group method with arithmetic mean. 

The system for searching for relationships between base sequences according to 
the present invention may further comprise an amplified sequence recognizing portion 
which reads pre-amplification base sequence data, and which produces post-amplification 
base sequence data based on primer data including recognition site data of primers used 
for the amplification, and the theoretical value calculating portion may calculate the 
theoretical restriction fragment length values based on the post-amplification base 
sequence data. Here, the term "amplification" means extraction of a certain portion of a 
base sequence located between a predetermined forward primer and a predetermined 
reverse primer. 

A method for searching for relationships between base sequences in genes 
according to the present invention comprises: a theoretical value calculating step for 
calculating theoretical restriction fragment length values of known genes digested by 
restriction enzymes based on restriction enzyme data in which kinds of restriction 
enzymes, recognition sites which are recognized by the restriction enzymes, and 
restriction sites at which the restriction enzymes cut in the recognition sites, are stored in a 
related manner, and base sequence data in which kmds of known genes and base 
sequences of the known genes are stored in a related manner, and outputting theoretical 
restriction fragment patterns in which kinds of known genes, the kinds of restriction 
enzymes, and theoretical restriction fragment length values posited to be obtained as a 
result of digesting known genes by the resfriction enzymes are stored in a related manner, 
and a comparing step for comparing the theoretical restriction fragment patterns with 
measured restriction fragment patterns including measured restriction fragment length 
values obtained as a resuh of digesting a sample by the restriction enzymes followed by 
measurement are stored in a related marmer, calculating the degree of similarity of the 
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known genes to the sample, and outputting it as analysis result data. 

The method for searching for relationships between base sequences according to 
the present invention may further comprise a displaying step for diagrammatically 
displaying samples having a high degree of similarity therebetween or a known gene and a 
sample having a high degree of similarity therebetween in a related manner based on the 
analysis result data. 

In the method for searching for relationships in base sequences according to the 
present invention, the displaying step may display samples having a high degree of 
similarity therebetween or a known gene and a sample having a high degree of similarity 
therebetween in a related manner by a dendrogram. 

In the method for searching for relationships between base sequences according 
to the present invention, the comparing portion may calculate the degree of similarity 
based on the theoretical restriction fragment patterns and the measured restriction 
fragment patterns by using an xmweighted-pair-group method with arithmetic mean. 

The method for searching for relationships between base sequences according to 
the present invention may further comprise an amplified sequence recognizing step for 
reading pre-amplification base sequence data, and producing post-amplification base 
sequence data based on primer data including recognition site data of the primer used for 
the amplification, and the theoretical value calculating step may calculate the theoretical 
restriction fragment length values based on the post-amplification base sequence data. 

A recording medium readable on a computer according to the present invention 
has recorded therein a computer program for implementing: a theoretical value calculating 
step for calculating theoretical restriction fragment length values of knovms gene digested 
by restriction enzymes based on restriction enzyme data in which kinds of restriction 
enzymes, recognition sites recognized by the restriction enzymes, and restriction sites at 
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which the restriction enzymes cut in the recognition sites are stored in a related manner, 
and base sequence data in which kinds of known genes and base sequences of the known 
genes are stored in a related maimer, and outputting theoretical restriction fragment 
patterns in which the kinds of known genes, the kinds of restriction enzymes, and 
theoretical restriction fragment length values are stored in a related manner; and a 
comparing step for comparing the theoretical restriction fragment patterns with measured 
restriction fragment patterns in which the kinds of restriction enzymes and measured 
restriction fragment length values obtained as a result of digesting a sample by the 
restriction enzymes followed by measurement are stored in a related maimer, calculating 
the degree of similarity of the known genes to the sample, and outputting it as analysis 
result data. 

The recording medium readable on a computer according to the present invention 
may further have recorded therein a computer program for implementing a displaying step 
for diagraimnatically displaying samples having a high degree of similarity therebetween 
or a known gene and a sample having a high degree of similarity therebetween in a related 
manner based on the analysis result data. 

In the recording medium readable on a computer according to the present 
invention, the displaying step may display samples having a high degree of similarity 
therebetween or a known gene and a sample having a high degree of similarity 
therebetween in a related manner by a dendrogram. 

In the recording medium readable on a computer according to the present 
invention, the comparing portion may calculate the degree of similarity based on the 
theoretical restriction fragment patterns and the measured restriction fragment patterns by 
using an unweighted-pair-group method with arithmetic mean. 

The recording medium readable on a computer according to the present invention 



7 



may further have recorded therein a computer program for implementing an ampHfied 
sequence recognizing step for reading pre-amplification base sequences, and producing 
post-amplification base sequences based on primer data including the recognition site data 
of primers used for the amplification, and the theoretical value calculating step may 
calculate the theoretical restriction fragment length values based on the post-amplification 
base sequence data. 

BRIEF DESCRIPTION OF DRAWINGS 
Fig. 1 is a block diagram showing the construction of the system for searching for 
relationships between base sequences according to a first embodiment of the present 
invention. 

Fig. 2 shows the first half of the post-amplification base sequence data according 
to the first embodiment. 

Fig. 3 shows the second half of the post-amplification base sequence data 
according to the first embodiment. 

Fig. 4 shows the structure of the restriction enzyme data according to the first 
embodiment. 

Fig. 5 shows an example of a user interface that orders calculation of a theoretical 
restriction fragment length value to the theoretical value calculating portion according to 
the first embodiment. 

Fig. 6 shows the structure of the theoretical restriction fragment patterns 
according to the first embodiment. 

Fig. 7 is a desktop image diagram showing a user interface for a measured value 
inputting portion according to the first embodiment. 

Fig. 8 shows the structure of the measured restriction fragment patterns according 
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to the first embodiment. 

Fig. 9 is a reference diagram representing a band distribution of fragments giving 
a basis for clustering a base sequence according to the first embodiment. 

Figs. lOA to IOC are tables showing the process of performing cluster analysis 
using a mean distance method based on the degree of similarity between base sequences 
according to the first embodiment of the present invention. 

Fig. 11 is a reference diagram showing the results of cluster analysis according to 
the first embodiment. 

Fig. 12 is a desktop image diagram showing an example according to the first 
embodiment where the displaying portion displays the results of analysis as a dendrogram. 

Fig. 13 is a block diagram showing the construction of the system for searching 
for relationships in base sequences according to a second embodiment of the present 
invention. 

Fig. 14 shows the first half of the pre-amplification base sequence data that serve 
as input data for the system for searching for relationships between base sequences 
according to the second embodiment. 

Fig. 15 shows the second half of the pre-amplification base sequence data that 
serve as input data for the system for searching for relationships between base sequences 
according to the second embodiment. 

Fig. 16 is a desktop image diagram showing an inputting picture for inputting a 
primer sequence for the recognition by an amplified sequence recognizing portion 
according to the second embodiment. 

Fig. 17 is a dendrogram showing an example of results obtained by performing 
classification and identification of species of bacteria using the system for searching for 
relationships between base sequences according to the present invention. 
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Figs. 18 to 27 are each a dendrogram showing an example of results obtained by 
performing classification and identification of species of bacteria using the system for 
searching for relationships between base sequences according to the present invention. 

Fig. 28 is a dendrogram showing combined results obtained by performmg 
classification and identification of species of bacteria using the system for searching for 
relationships in base sequences according to the present invention. 

Fig. 29 is a table showing results of the search for a species name of an isolated 
denitrifying bacterium by the method of the present invention in which a species name of 
the isolated bacterium was posited, with respect to a representative strain the partial 
sequence of its 16S rDNA was determined, and a homology was searched v^th reference 
to base sequences in the open database to find the species name of the bacterium having 
the closest base sequence structure of 16S rDNA. 

DETAILED DESCRIPTION OF THE INVENTION 
Hereinafter, embodiments according to the present invention will be explained in 
detail with reference to the attached drawings. 
First Embodiment 

Fig. 1 is a block diagram showing the construction of the system for searching for 
relationships between base sequences according to a first embodiment of the present 
invention. In Fig. 1 , restriction enzyme data 2 1 shows information on restriction 
enzymes appropriately set, which includes kinds of restriction enzymes, the recognition 
sites which are recognized by the restriction enzymes, and restriction sites at which the 
restriction enzymes cut in the recognition sites, in a related manner. The restriction 
enzyme data 21 may contain data on one or more restriction enzymes. 

A post-amplification base sequence data 22b shows information on a knovra gene 
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in a region where it is expected to be amplified by a polymerase chain reaction using the 
known gene as a template and using a desired primer set (including a forward primer and a 
reverse primer). The post-amplification base sequence data 22b includes the kind of 
known gene and the base sequence of the known gene in the region to be amplified in a 
related manner. The number of the above-described primer set and known gene may be 
one or more and the above-described post-amplification base sequence data 22b may show 
information on one or more known genes in one or more kinds of regions. The 
post-amplification base sequence data 22b can be obtained fi-om the open base sequence 
database, and be manually or otherwise edited, for example. 

In a theoretical value calculating portion 11, theoretical restriction Augment 
length values of known genes, which are expected to be obtained by digesting the known 
genes in the amplified regions by the restriction enzymes, are calculated based on the 
above-described restriction enzyme data 21 and post-amplification base sequence data 22b, 
and theoretical restriction Segment patterns 23 are output. 

hi a measured value inputting portion 12, measured restriction fi-agment patterns 
25 are input. The measured restriction fi-agment patterns 25 include measured restriction 
fi-agment length values obtained by digesting a sample by the same restriction enzymes as 
used for calculating the above-described theoretical restriction fi-agment lengtii values and 
measuring the lengths of the digested base sequences as well as the kind of the restriction 
enz5Tiies in a related manner. 

In a comparing portion 13, the theoretical restriction fragment patterns 23 and the 
measured restriction firagment patterns 25 are compared with each other, the degree of 
similarity between the theoretical restriction fragment patterns 23 and the measured 
restriction fi-agment patterns 25 is calculated, and the degree of similarity of the known 
gene and the sample is output as analysis result data 27. 
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An allowance setting portion 14 sets up allowance limits of error for comparing 
the theoretical restriction fragment patterns 23 and the measured restriction fragment 
patterns 25. The allowance limits of error are stored as allowance limits data 29. The 
comparing portion 13 refers to the allowance limits data 29 and treats the theoretical 
restriction fragment patterns 23 and the measured restriction fragment patterns 25 as being 
identical with each other when the difference therebetween is within the allowance limits 
of error, and then calculates the degree of similarity between the theoretical restriction 
fragment patterns and the measured restriction fragment patterns 25. 

A displaying portion 15 relates samples having a high degree of similarity 
therebetween to each other or relates a known gene and a sample having a high degree of 
similarity therebetween to each other based on the analysis result data 27 output by the 
comparing portion 13 and displays them diagrammatically. 

Here, the details of the post-amplification base sequence data 22b will be 
described. Fig. 2 shows the first half of the post-amplification base sequence data 22b 
(base sequence data) while Fig. 3 shows the second half of the post-amplification base 
sequence data 22b. As shown in Figs. 2 and 3, the post-amplification base sequence data 
22b are text format data. The first half of the post-amplification base sequence data 22b 
includes at least the DNA number (ACCESSION), the name of organism of origin 
(ORGANISM), and the name of the known gene (KEYWORDS) while the second half of 
the post-amplification base sequence data 22b includes known gene base sequence 
(ORIGIN) (SEQ ID NO. 1) in the region where it is expected to be amplified by 
polymerase chain reaction using the known gene as a template and using a desired primer 
set (forward primer and reverse primer). In the base sequence data, the letter "a" 
represents adenine, "g" represents guanine, "c" represents cytosine, and "t" represents 
thymine. 
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Next, the details of the restriction enzyme data 21 will be described. Fig. 4 
shows the data structiire of the restriction enzyme data 2 1 . The restriction enzyme data 
2 1 are registered in the restriction enzyme databeise, which are shown in a table format 
having columns for the restriction enzyme nimiber, the restriction enzyme name, the 
recognition site recognized by each restriction enzyme, and the restriction site at which the 
restriction enzyme cut in the recognition site, as shown in Fig. 4. The restriction enzyme 
number is a number automatically given when a new restriction enzyme data is added to 
the database through a user interface (not shown in the figures). The recognition site 
recognized by each restriction enzyme includes the restriction site at which the restriction 
enzyme cut in the recognition site. The restriction site is indicated by the number of 
bases from the 5'-terminal side of the recognition site. For example, Smal, which is set 
to be restriction enzyme number 12 in Fig. 4, recognizes a recognition site . .cccggg. . . " 
of a target DNA molecule, and cuts it at the restriction site which is the third bond, 
between the third base (c) and the fourth base (g). As a result, the target DNA molecule 
is digested into "... ccc" and "ggg. . . ". At present, data on 57 kinds of restriction enzymes 
are registered. 

Fig. 5 is a desktop image diagram of an example of a user interface that orders 
the calculation of the theoretical restriction firagment length value in the theoretical value 
calculating portion 11. In a central portion of the user interface, there are displayed lists 
showing kinds of known genes contained in the post-amplification base sequence data 22b 
and the restriction enzyme names contained in the restriction enzyme data 21, respectively. 
On the lower portion of the user interface are displayed a "RESTRICTION FRAGMENT 
LENGTH CALCULATION" button and a "RETURN" button. 

The user may select the desired kind of known gene and restriction enzyme name 
and push the "RESTRICTION FRAGMENT LENGTH CALCULATION" button on the 
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lower portion of the user interface to give the instruction for the calculation of the 
theoretical restriction fragment length value to the theoretical value calculating portion 1 1 . 

When the "RESTRICTION FRAGMENT LENGTH CALCULATION" button is 
pushed, the theoretical value calculating portion 1 1 reads the base sequence of a known 
gene in the amplified region of the post-amplification base sequence data 22b and the 
recognition site of the restriction en2yme from the restriction enzyme data 21, and 
compares them with each other. A search is conducted from the 5' end terminal of the 
base sequence of the known gene to determine whether or not the recognition site exists in 
the base sequence of the known gene in the region to be amplified. If recognition sites 
are present in the base sequence of the known gene, the number of bases present from the 
recognition site closest to the 5' end terminal of the base sequence up to the 5' end terminal 
is calculated as a theoretical resfriction fragment length value. Then, the numbers of 
bases present between the adjacent restriction sites are respectively calculated and are 
input as theoretical restriction fragment length values. Then, the theoretical restriction 
fragment length values are output as theoretical restriction fragment pattern 23. 

Fig. 6 is a table showing an example of the structure of theoretical restriction 
fragment patterns 23. As shown in Fig. 6, the theoretical restriction fragment patterns 23 
are registered in the restriction fragment database, which are shown in a table format 
having columns for indicating the items of DNA number, numerical value section, name 
of organism of origin, name of known gene, restriction enzyme number, and restriction 
fragment length. Further, the theoretical restriction fragment pattern 23 may have further 
a column indicating data for distinguishing digested fragments (base sequence from the 
restriction site closest to the 5' end terminal of the base sequence up to the 5' end terminal, 
or base sequence between the adjacent restriction sites). It is possible to distinguish 
digested fragments having the same restriction fragment length obtained by using the 
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same restriction enzyme from each other. The theoretical restriction fragment patterns 23 
are tabularized for every DNA (distinguished by DNA number), for every restriction 
enzyme (distinguished by restriction enzyme number) and for each restriction fragment 
length value. The DNA numbers, the names of organisms of origin, and the names of 
known genes are obtained from the post-amplification base sequence data 22b and are 
output by the theoretical value calculating portion 1 1 . The restriction enzyme numbers 
are read from the restriction enzyme data 21 and are output by the theoretical value 
calculating portion 1 1 . The theoretical restriction fragment length values are calculated 
by the theoretical value calculating portion 1 1 by the above-described procedures. The 
numerical value section indicates the classification of theoretical values or measured 
values. The numerical section " 1 " indicates that the restriction fragment length is a 
theoretical value. 

Next, the measured restriction fragment pattern 25 will be described. 

To measure the restriction fragment length of a sample, the base sequence 
contained in the sample is amplified by polymerase chain reaction (PGR). As a primer 
set used for the PGR, the primer set used for obtaining the post-amplification base 
sequence data 22b is preferred. After the amplification, the amplified DNA molecule is 
digested by the restriction enzyme used for obtainmg the theoretical restriction fragment 
pattern 23 . Then, the lengths of digested fragments are measured, for example, by an 
agarose gel elecfrophoresis method, polyacrylamide gel electrophoresis method, a 
microchip electrophoresis method, or the like method, and are input into the measured 
value inputting portion 12 as measured restriction fragment length value. 

Fig. 7 is a desktop image diagram showing an example of a user interface of the 
measured value inputting portion 1 2. In the user interface shown in Fig. 7, the user 
inputs the DNA number, selects the name of a restriction enzyme from the list of 
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restriction enzyme names displayed based on the restriction enzyme data 21, and records 
the measured restriction fragment length value so that the measured restriction fragment 
patterns 25 can be input into the measwed value inputting portion 12. Furthermore, the 
user can group a plurality of samples, give group names to the resulting groups, and then 
input the group names into the measured value inputting portion 12. When the microchip 
electrophoresis device is used in order to obtain the above-described measured restriction 
fragment length values, the system may be arranged such that the measured restriction 
fragment length values can be automatically input into the measured value inputting 
portion 12 from the microchip electrophoresis device. 

Fig. 8 is a table showing an example of the data structure of measured restriction 
fragment patterns 25. The measured restriction fragment patterns 25 are registered in the 
measured restriction fragment database, which are shown in table format in Fig. 8 having 
columns indicating the items of the DNA number, the numerical value section, the name 
of the organism, the name of the gene, the restriction enzyme number, and the measured 
restriction fr^ment length value. The measured restriction fragment patterns 25 are 
tabularized for every DNA (distinguished by DNA number), for every restriction enzyme 
(distinguished by restriction enzyme number), and for each measured restriction fragment 
length value. The numerical value section indicates the classification of theoretical 
values or measured values. The numerical section "2" indicates that the restriction 
fragment length is a measured value. 

The comparing portion 13 compares the theoretical restriction fragment patterns 
23 with the measured restriction fragment patterns 25, and analyzes the relationship 
therebetween. In the present embodiment, the comparing portion 13 produces analysis 
result data 27 by using an unweighted-pair-group method with arithmetic mean 
(UPGMA). 
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Fig. 9 is a reference diagram representing the measured restriction fragment 
pattern or the theoretical restriction fragment pattern, which gives a basis for clustering 
base sequences by the comparing portion 13. Fig. 9 is a simplified image of the 
restriction pattern obtained by electrophoresis. In Fig. 9, symbols A, B, and C each 
represent respective DNA numbers contained in the theoretical restriction fragment 
patterns 23 or measured restriction fragment patterns 25. The bands in 3 rows m the 
longitudinal direction indicate patterns of restriction fragments obtained from known 
genes or samples distinguished by the DNA numbers by bands arranged in accordance 
with the theoretical restriction fragment length value or measured restriction fragment 
length values in a descending order downwardly. Also, in Fig. 9, numerals from 1 to 8 
indicate band numbers for distinguishing the respective bands. 

The comparing portion 13 reads the theoretical restriction fragment patterns 23 or 
the measured restriction fragment patterns 25 and calculates the degree of similarity 
between the data given by two DNA numbers. For example, the definition of degree of 
similarity S(A, B) is given by the following formula: 

S(A, B) (total number of common bands that A and B have)x2/(number of bands that A 
has + number of bands that B has). 

In Fig. 9, A has five bands with band numbers of 1, 2, 4, 5 and 8, respectively, 
while B has six bands with band numbers of 1 , 3, 4, 5, 7, and 8, respectively. Therefore, 
the four bands with band numbers 1, 4, 5, and 8 are common bands that A and B have, and 
the total number of common bands that A and B have is 8. Therefore, 

S(A, B) = 4x2/(6+5) 
and hence the degree of similarity is calculated to be 8/11 (about 0.727). 

As the measured restriction fragment length value contains an error, the 
comparing portion 13 uses allowance limits data 29 set by an allowance setting portion 
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before a judgment can be made as to whether or not a certain band is a common band. 
More particularly, when the distance betvveen the bands to be judged is within the 
allowance limits of error contained in the allowance limits data 29, they are treated as 
common bands. In Fig. 9, the allowance limits of error is expressed as a ratio of the 
longitudinal distance between the bands compared with each other to the longitudinal 
distance between the bands with band numbers 1 and 8. The allowance limits of error is 
properly set depending on the electrophoresis method used and preferably is on the order 
of 10% to 5%. For example, preferably it is about 10% in agarose gel electrophoresis, 
about 7% in low melting point agarose gel electrophoresis, and about 5% in 
polyacrylamide gel electrophoresis. If the allowance limits of error is set too low, a 
problem arises in that inherently common bands will appear to be different bands. On 
the other hand, if the allowance limits is too high, a problem arises in that inherently 
different bands will appear to be common bands. 

Next, the comparing portion 13 performs cluster analysis of base sequences based 
on the calculated degrees of similarity. Figs. 1 OA, 1 OB, and 1 OC are tables showing an 
example of the process of cluster analysis performed by the comparing portion 13 using a 
mean distance method based on the degrees of similarity. In Figs. 1 OA, 1 OB, and 1 OC, 
four base sequences A, B, C, and D are taken as subjects of analysis. Fig. lOA shows the 
calculated degrees of similarity between any two of A, B, C, and D. In the table shown 
in Fig. lOA, the combination that shows the highest degree of similarity is A and C (the 
degree of similarity: 0.8). Accordingly, A and C are combined to form a cluster of 
(A+C) and the next step is followed. 

Fig. lOB shows the degrees of similarity calculated by substituting A and C in 
Fig. 1 OA by a single cluster of (A+C). 

For example, as the similarity between B and the cluster (A+C), i.e., S(B, A+C), 
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a mean value of S(B, A) and S(B, C) is calculated to be 0.6. Also, as the degree of 
similarity between D and the cluster (A+C), i.e., S(D, A+C), a mean value of S(D, A) and 
S(D, C) is calculated to be 0.2. Then, in the table shown in Fig. 1 OB, the combination 
that shows the highest degree of similarity is (A+C) and B. Accordingly, (A+C) and B 
are combined to form a cluster of ((A+C)+B), and the next step is followed. 

Fig. IOC shows the degrees of similarity calculated by substituting (A+C) and B 
in Fig. 1 OB by a single cluster of ((A+C)+B). Here, as the similarity between ((A+C)+B) 
and B, i.e., S(D, (A+C)+B), a mean value of S(D, A+C) and S(D, B) is calculated to be 
0.3. 

Fig. 1 1 is a reference diagram showing the analysis result data obtained by the 
unweighted-pair-group method with arithmetic mean as explained above in a dendrogram. 
In Fig. 11, the numerical value attached to the node portion of the tree represents the 
degree of similarity between the terminal nodes (A, C, B, or D) or the non-terminal nodes 
(cluster (A+C), ((A+C)+B)). For example, the degree of similarity between B and the 
cluster (A+C) is calculated to be 0.6. 

The displaying portion 15 reads the analysis result data 27 and diagrammatically 
displays the samples registered by the user and a known gene having a high degree of 
similarity thereto in a related manner. Fig. 12 is a desktop image diagram showing an 
example of an analysis result displayed by the displaying portion 15 as a dendrogram. 
Fig. 12 shows the results of cluster analysis of a plurality of DNA samples (API, AP2, and 
APS) and a plurality of known genes. 

Second Embodiment 

In the second embodiment, the post-amplification base sequence data 22b are 
automatically produced based on the pre-amplification base sequence data 22a. In other 
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words, the region that is to be amplified by the polymerase chain reaction using a primer 
set (forward primer and reverse primer) is automatically recognized in the 
pre-amplification base sequence of a known gene and the post-amplification base 
sequence is extracted. 

Fig. 13 is a block diagram showing the construction of the system for searching 
for relationships in base sequences according to the second embodiment. In Fig. 13, the 
pre-amplification base sequence data 22a include the infonnation of the pre-amplification 
base sequence of known genes. An amplified sequence recognizing portion 3 1 reads the 
pre-ampUfication base sequence data 22a, recognizes the base sequence of the region to be 
amplified by using the primer set based on the primer data including the base sequence 
information of the primer set, and outputs the base sequence of the region as the 
post-amplification base sequence data 22b. 

Fig. 14 shows the first half of the pre-amplification base sequence data 22a while 
Fig. 15 shows the second half of the pre-amplification base sequence data 22a. In Fig. 15, 
the second half of the pre-amplification base sequence data 22a includes the base sequence 
information of a known gene before amplification (SEQ ID NO. 2). As shown in Fig. 15, 
the pre-amplification base sequence data 22a are text format data similar to the 
post-amplification base sequence data 22b. Fig. 3 and Fig. 1 5 relate to the same gene 
(ACCESSION="M59070", KEYWORDS=" 1 6S ribosomal RNA.") The 
post-amplification base sequence of the known gene shown in Fig. 3 (SEQ ID NO. 1) 
corresponds to the 16S rDNA sequence from 23 to 1031, which represents the 
pre-amplification base sequence of the 16S rDNA gene shown in Fig. 15 (SEQ ID No. 2) 
and its length is 1009 bp. 

At fnst, the primer data, which includes information on the base sequences, the 
names, and the allowance limits for mismatch of predetermined forward primer and 
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reverse primer, are input and registered in the primer database. The folder having the 
names of both primers, in which the post-amplification base sequence is saved, is 
automatically formed when a new primer data is added to the database (not shown in the 
figures). Fig. 16 is a desktop image diagram showing a selection of the primer data of 
the primer set as described above. (In the screen shown in Fig. 16, the user can select the 
respective primer data.) The forward primer used in Fig. 16 is named "41f ', has a base 
sequence of "gctcagattgaactcggcg" of SEQ ID NO. 3, and has an allowance limit for 
mismatch of 4%. Similarly, the reverse primer is named "1066r", has a base sequence of 
"acatttcacaacacgagctg", and has an allowance limit for mismatch of 4. The primer data 
thus selected are utilized by the amplified sequence recognizing portion 3 1 when the base 
sequence in the region to be amplified by using the primer set is extracted from the 
pre-amplification base sequence of the known gene. 

The amplified sequence recognizing portion 31 reads the pre-amplification base 
sequence data 22a, and compares the pre-amplification base sequence of the known gene 
with the base sequence of the forward primer. Within the allowance limits for mismatch, 
the region, where the pre-amplification base sequence and the base sequence of the 
forward primer match with each other, is searched from the 5' terminal side thereof in 
order. Then, the pre-amplification base sequence is compared with the base sequence of 
the reverse primer. Within the allowance limits for mismatch, the region, where the 
pre-amplification base sequence and the base sequence of the reverse primer match with 
each other, is searched fi-om the 3' terminal side thereof in order. The base sequence in 
the region sandwiched by the forward primer and the reverse primer is output as the 
post-amplification base sequence of the known gene (ORIGIN) in the post-amplification 
base sequence data 22b and is saved in the primer folder described above. 

The above-described system for searching for relationships between base 
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sequences in genes is realized by using a computer system. A computer program that 
enables a computer to implement respective operations of the above-described theoretical 
value calciilating portion, measxured value inputting portion, comparing portion, allowance 
limits setting portion, and displaying portion is stored in a recording medium readable on a 
computer. The computer reads out the computer program and implements it so that the 
above-described operations can be performed. Here, the recording medium readable on a 
computer includes a floppy disk, a magneto-optical disk, CD-ROM, DVD-ROM, a 
magnetic hard disk, a semiconductor memory and the like. 

As the computer system for realizing the system for searching for relationships 
between base sequences in genes, a general-purpose computer such as a personal 
computer or a workstation can be used. 

The format of data treated by the system for searching for relationships between 
base sequences in genes is not limited to those described in the first and second 
embodiments described above, and data in other formats or expressions may also be used. 
A portion of the data in a table format as described in the above embodiments is in a 
non-normalized form. Such data may be normalized. For example, in the theoretical 
restriction fragment patterns shown m Fig. 6, the DNA numbers, the names of the 
organisms, and the names of genes are shown in the same table. However, they may be 
shown in separate tables. 

Next, test examples in which classification of bacteria using the system for 
searching for relationships between base sequences in genes according to the present 
invention are performed will be described. 

Test Example 1 

The 16S rDNA sequence data is downloaded genus by genus from the DNA 
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information supplied by the National Institute of Agrobiological Sciences, Ministry of 
Agriculture and Fisheries, Japan. The DNA information supplied by the National 
Institute of Agrobiological Sciences is constructed based on the DNA information 
obtained from the Gene Bank (National Center for Biotechnology Information), DDBJ 
(National Institute of Genetics), and EMBL (European Molecular Biology Laboratory). 
Also, 5'-gctcagattgaacgctggcg-3' as the forward primer (41f), 5'-acatttcacaacacgagctg-3' as 
the reverse primer (1066r), and fourteen kinds of restriction enzymes were input. The 
region sandwiched by the forward primer and the reverse primer was extracted as 
post-amplification base sequence data 22b from the pre-amplification base sequence data 
22a by the amplified sequence recognizing portion 31, and 357 genera, 1233 species, or 
1503 kinds of base sequence data were obtained as the post-amplification base sequence 
data 22b. Then, the lengths of fragments digested by the restriction enzymes were 
calculated and output as theoretical restriction fragment patterns by the theoretical value 
calculating portion 1 1 , 

Then, alkylphenol decomposing bacteria isolated from various soils (from 
Kyushu Kyoritsu University, Meijo University) were used as samples API to AP9 and 
alkali resistant bacteria (from Yamaguchi University) were used as samples MAI to MA4. 
Respective chromosomal DNAs was extracted from the samples by a conventional method 
and PGR was performed by using the forward primer and reverse primer to amplify 16S 
rDNA. The amplified 16S rDNA was digested by the restriction enzymes used for 
calculating the theoretical restriction fragment length value. Then, the lengths of the 
resulting digested fragments were measured using densitometory and defined as measured 
restriction fragment length values. As sample MAU, the same strain as the sample MAI 
was used, the same operations as for MAI were performed, the measured restriction 
fragment length values were determined, and measured restriction fragment patterns were 
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obtained. 

Fig. 17 is a dendrogram showing the relationship between the base sequences 
obtained from each sample output as analysis result data 27 by the comparing portion 13 
based on the measxired restriction fragment patterns obtained as described above. 
Also, each of Figs. 18 to 27 is a dendrogram showing the relationship between the base 
sequence obtained from each sample and the base sequence of a known gene output as the 
analysis result data 27 by the comparing portion 1 3 . Figs. 1 8 , 1 9, 20, 2 1 , 22, 23 , 24, 25 , 
26, and 27 show the results of API, AP3, AP6, AP2, AP5, MAI, MAll, MA2, MA3, and 
MA4, respectively. In analysis resuh data 27 by the comparing portion 13, the 5 known 
sequences, which have the restriction fragment patterns of the highest similarity with the 
each measured restriction fragment patterns, are automatically searched in theoretical 
restriction fragment patterns 23 (total 1503 sequence), and shown in the Figs. 18 to 27. 

From the above-described analysis results, it was revealed that samples MAI and 
MAll have the closest relationship to Aeromonas hydrophila of the Aeromonas group 
belonging to the Proteobacteria c subdivision and also that they have close relationships to 
the genus Serratia, the genus Shewanella, or the genus Plesiomonas of the 
Enterobacteriaceae belonging to the Proteobacteria c subdivision. 

Also, it was revealed that the sample MA2 even has a close relationship to the 
genus Sporolactobacillus or the genus Bacillus (both are low GC gram positive bacteria, 
Bacillus/Clostridium, and Bacilluceae (Bacillus family)). 

The sample MA3 was revealed to have the closest relationship to Bacillus cereus 
and Bacillus thuringiensis (low GC gram positive bacteria, Bacillus/Clostridium, and 
Bacilluceae). Note that Bacillus cereus and Bacillus thuringiensis have been indicated to 
be related to each other and there are many opinions that they are of the same species. 
The sample MA4 was indicated to have the closest relationship to the genus 
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Oeskovia or the genus Cellulomonas (both are Actinobacteria, Actinobacteriadae, 
Actinomycetales, Microccccineae, and Cellulomonadaceae) and a close relationship to the 
genus Actinomycetes (Actinobacteria, Actinobacteridae, Actinomycetales, 
Actinomyccineae, and Actinomycetaceae), and the genus Streptomyces (Actinobacteria, 
Actinobacteridae, Actinomycetales, Streptomycineae, and Streptomycetaceae). 

It was revealed that the samples API and AP3 each have a close relationship to 
Pseudomonas putida, Pseudomonas fulva, Pseudomonas straminea, Pseudomonas 
alcaligenes, and Flavimonas oryzihabitans (each belonging to the Proteobacteria c 
subdivision, Pseudomonadaceae). 

Further, the sample AP6 has a close relationship to bacteria of the genus 
Xanthomonas (Proteobacteria c subdivision, Lysobacterrales, and Xanthomonas group) 
and the samples AP2 and APS were posited to belong to the genus Phyllobacterium, the 
genus Rhizobium, the genus Agrobacterium (Proteoacteria, a subdivision, Rhizobiaceae 
group, and Phylobacteriaceae or Rhizobiaceae, respectively). 

Fig. 28 is a dendrogram showing the relationship between the base sequence 
obtained from the entire sample and the base sequence of the known gene output as the 
analysis result data 27 by the comparing portion 13. Analysis results of each sample 
shown in Figs. 18 to 27 can be automatically combined, which enable analysis of many 
samples simultaneously in a short time. 

Thus, using the present test example, data on very many kinds of known genes 
including base sequence data newly included in the futvire could be obtained at low cost 
and in a short time based on publically available data. Upon performing classification 
and identification of the kind of unknown gene contained in the sample, comparison with 
many kinds of known genes can be performed at high speed so that the degree of 
similarity between known bacterial strains and unknown bacterial strains can be calculated 
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efficiently and with high precision. 
Test Example 2 

The species names of 134 strains of denitrifying bacterium isolated by a 
conventional method were inferred based on the above-described method. Thereafter, 
the partial 16S rDNA sequences of representative bacterial strains were actually 
determined by the ABI 373 sequencing system (Perkin-Elmer Japan). Then, a homology 
search was performed using public sequence databases using the BLAST algorithm, so as 
to search for the species name of the bacterium having the closest base sequence to that of 
the 16S rDNA. Fig. 29 shows the group nimber (I to XIII) obtained by classifying the 
134 strains of isolated denitrifying bacterium, the number of isolates contained in the 
group, the species name inferred by the above-described method (RFLP), the name of a 
bacterium having a base sequence having a high degree of homology revealed by actually 
determining the base sequence of 1 6S rDNA, and degrees of difference. The degrees of 
difference as described above were obtained by determining the 16S rDNA sequence 
obtained from each isolate belonging to the same group, performing a homology search 
for the obtained base sequences, respectively, indicating the ratio of the number of 
matched bases to the total number of bases in the base sequence in the region having high 
homology in terms of percentage (%). This revealed that the species names posited by 
the above-described method well correspond with those obtained by the comparative 
analyses of base sequences. 

As described above, according to the present invention, it is possible to calculate 
the degree of similarity between a known gene and a sample efficiently and with high 
precision. Fiorthermore, use of the present invention enables not only a search of the 
relationships between procaryotic organisms, but also a search of the relationships 
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between eucaryotic organisms. Also, in order to relate those eucaryotic organisms 
having a high degree of similarity to each other based on the degree of similarity between 
many kinds of known genes and a sample, and diagrammatically displaying the 
relationship using a dendrogram, it is possible to provide analysis results in a form that is 
understandable to users. 

Furthermore, according to the present invention, the amplified sequence 
recognizing portion automatically recognizes the amplified portion and produces 
post-amplification base sequence data and hence post-amplification base sequence data 
can be prepared without much time and labor so that the analysis of base sequences can be 
made more efficient. 

Moreover, according to the present invention, realization of the system for 
searching for relationships in base sequences using a general purpose computer such as a 
personal computer or a workstation can provide a system that is extremely inexpensive 
compared with a special sequencer or the like. 



